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Abstract 

QCD at the classical level possesses scale invariance which is broken by quantum effects. 
This "dimensional transmutation" phenomenon can be mathematically described by formu- 
lating classical gluodynamics in a curved, conformally flat, space-time with non-vanishing 
cosmological constant. We study QCD high-energy scattering in this theory. We find that 
the properties of the scattering amplitude at small momentum transfer are determined by 
the energy density of vacuum fluctuations. The approach gives rise to the power growth of 
the total hadron-hadron cross section with energy, i.e. the Pomeron. The intercept of the 
Pomeron and the multiplicity of produced particles are evaluated. We also speculate about 
a possible link between conformal anomaly and TeV scale quantum gravity. 
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1 Introduction 



At the classical level, QCD is a scale invariant theory. Quantum fluctuations, however, give rise 
to scale anomaly Q which is associated with dimensionful scale A ~ 0.2 GeV. The perturbative 
expansion is possible if the typical distance Q" 1 at which the interaction occurs is much smaller 
than A -1 , so the coupling is weak @: a s (Q) <C 1. However, this does not mean that at short 
distances non-perturbative contributions can always be neglected. Indeed, the strength of the 
semi-classical vacuum fluctuations of the gluon field is inversely proportional to the coupling 
constant and thus can be large at short distances, when ~ l/a s (Q) S> 1. The goal of the 
present work is to derive an effective lagrangian for gluodynamics which makes it possible to 
systematically treat such contributions to the scattering amplitudes. 
Gluodynamics is defined by the lagrangian 

£gluon = ~~^F f±v F flv ■ (1) 

If treated classically, this theory is invariant under the scale transformations x —* Ax. Transfor- 
mation of scale is associated with the scale current s^; its divergence is equal to the trace of the 
energy-momentum tensor: 

a% = o*(x). (2) 

Without quantum effects, 0^{x) = 0, and the theory is scale invariant. Quantum effects break 
scale invariance However, the broken symmetry still manifests itself in the following set of 
low energy theorems for different Green functions involving operator 0fi(x) which can be proven 
by using renormalization group arguments ||: 

i n J dx 1 ...dx n (0\T{6^(x 1 ),...,9^{x n )}\0) connect ^ = <^(x)) vac (-4r (3) 

In Ref. Q it was shown that these low-energy theorems entirely determine the form of 
the effective low energy lagrangian for gluodynamics, which describes semi-classical vacuum 
fluctuations of gluon field at large distances. It is formulated in terms of an effective scalar 
dilaton field: 

/:dilaton = ^4^/ 2 (^ X ) 2 + kv|e X (l-x) , (4) 

where the real scalar field \ of mass m represents the dilaton, and — |e v | is the vacuum energy 
density. 

The idea underlying the derivation can be formulated in the way which illustrates the analogy 
between gluodynamics and general relativity: QCD coupled to the conformally flat gravity is 
scale and conformally invariant in any number of dimensions, and thus such theory is anomaly- 
free. However, the vacuum energy density of flat-space gluodynamics — |e v | (which plays a role 
of cosmological constant in general relativity) manifestly breaks scale and conformal symmetry 
of gluodynamics in conformally flat space-time. The low energy effective lagrangian (H) thus 
can be derived from the Einstein-Hilbert lagrangian for the described theory. Note that in 
the limit |e v | — ► lagrangian (|j) is scale and conformally invariant, since another dimensionful 
parameter m just sets the normalization of the dilaton field x- At small dilaton momenta 
p -C m the lagrangian (||) satisfies the low energy theorems (||). It is remarkable that the 
effective lagrangian involves only one scalar field. The reason is that the gravity is excited by 
the trace of energy momentum-tensor, so there is only one independent Einstein equation. 

In the present work we are going to extend the range of validity of the effective lagrangian 
(|j) to shorter distances by including the interactions of dilaton field with gluons. Generally, it 
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can be obtained from the lagrangian of gluodynamics by integrating out soft gluon modes up to 
the scale m which thus becomes the dimensionful scale of the theory. This implies that the scale 
A and the dilaton mass m are intimately connected. Actually, as we will show, it appears more 
convenient to construct the generalized effective lagrangian using Einstein equations once again. 
In the next section we employ this method to obtain the lagrangian (||). One can easily check 
that it is non-renormalizable. Therefore, we need to specify the additional dimensionful scale - 
the ultra-violet cutoff Mq. Obviously, Mq is a function of |e v | and m. It can be calculated using 
the requirement of vacuum stability, i.e. all radiative corrections must cancel out to give correct 
value of vacuum energy density. Its numerical value turns out to be quite large, Mq ~ 3 GeV. 
This is in agreement with Ref. ||, || ^ where it has been argued that the typical scale of the 
vacuum gluon density due to semi-classical fluctuations can indeed be that large. In sec. 2 we 
discuss the effective lagrangian (||) in detail and calculate the scale Mq. 

In sec. 3 we apply the effective lagrangian (||) to calculate a contribution of the non- 
perturbative semi-classical fluctuations of the QCD vacuum to the pomeron intercept. Ex- 
perimental data on different reactions at high energies and small momentum transfer support 
existence of the pole in the scattering amplitude responsible for high energy behavior of the total 
cross sections 0. This leads to the power-like behavior of the total cross section cr t ot = <is ap ~ l , 
where ? depends upon particular reaction and ap is the universal parameter called the inter- 
cept. After decades of hard work, understanding the nature of the object exchanged in t-channel 
(pomeron) is still a challenging problem. It turned out that the perturbative QCD is far from 
giving a clear answer even if the scales inherent to colliding particles are hard. Partly, this is 
due to diffusion of the transverse momenta in the partonic cascade towards the non-perturbative 
scales. This implies that the main contribution to the pomeron stems, perhaps, from the non- 
perturbative sector of QCD. Assuming that some non-perturbative fluctuations of vacuum dom- 
inate the scattering amplitude at high energies one has to explain the power-like behavior of 
the total cross section cr to t ~ s ap ~ l and to obtain the correct numerical value of ap. A few 
authors have reported on this problem. Contribution of instantons to the scattering amplitude 
at high energies was investigated in Ref. p. In Ref. [10] it has been argued that fluctuations 
of strings in the framework of AdS/CFT theory can yield the desired result. In Refs. [11, 12] 
the spectrum of glueballs in QCD was discussed as well as its relation to the structure of the 
pomeron. Possible contributions of extra compact dimensions was discussed in Ref. [13, 14]. A 
more phenomenological approach has been developed in framework of the BFKL gluon emission 
with additional assumption for such an emission in the non-perturbative QCD domain [15]. In 
Ref. jn| two of us have used the fact that the non-perturbative QCD vacuum is dominated by 
the semi-classical fluctuations of the gluon field which is inversely proportional to the coupling 
F ~ 1/g. The importance of such fluctuations can be understood in the following way. Consider 
a gluon ladder in each rung of which two gluons in a singlet state are produced in s-channel. 
Each rung gives contribution of the order g 2 F 2 ~ 0(1) since the two- gluon singlet state gives 
contribution proportional to the trace of energy-momentum tensor. Such ladder diagrams can 
be summed up yielding a reasonable value of the pomeron intercept. In the present paper we 
suggest a method for systematic calculation of processes involving such non-perturbative modes 
of gluodynamics. In particular, we will calculate the total cross section at high energies, repro- 
duce its power- like behavior and numerical value of the intercept. We also calculate the final 
state multiplicity in this approach. 

We conclude the paper by presenting some speculative ideas on the possible link between 
QCD and gravity. Recall that the effective lagrangian (|5|) can be derived from the General 
Relativity. Thus the gravitational constant can be expressed in terms of the dilaton mass and 
the vacuum energy density, Eq. di~l|). Its numerical value G = 0.2 GeV" 2 , of course, has nothing 
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to do with the Planck gravitational constant. Nevertheless, it may be possible to interpret this 
result as a signature of the strong gravitational interactions at scale much lower than the Planck 
one. If we accept this interpretation, then it turns out that the gravitational effects are essential 
already at the GeV scale, in order to be responsible for the QCD vacuum structure. Several 
authors have suggested mechanisms by which gravitational interactions can become strong at 
scales much larger than the Planck one [17|. Recently it has been suggested that quantum 



gravity can become strong at the TeV scale [18]. In sec. 4 we briefly discuss those ideas. 



2 Dilaton effective lagrangian 

The dilaton effective lagrangian in gluodynamics describing the propagation and interactions of 
gluon and dilaton fields can be constructed as a generalization of Eqs. (Q),(|l|). It reads 

£= Nl e x/2 (9 )2 + |ev|e x (1 _ x) _ e X (1 _ x) V F a^ (5) 

m z 2 4 p 

Eq. (|5|) can be obtained utilizing Legendre transformation of the lagrangian describing coupling 
of gluodynamics to gravity 

1 1 

8^G R ~ 4' 



C '= [^ R --aKu F<1 ^ + \^\) V- det ^ ■ (6) 



The first term in (jfj) corresponds to the gravitational lagrangian Cq- It is expressed through the 
Ricci scalar R which can be easily calculated for the conformally trivial metric g^ u = ij^e 11 ^ 

1 . u 1,„ ,„\ 13 



Cg = toG*N-**9,~ = ^ ^ h+ -( dtlh y) = -g^cw , (7) 

where we have used integration by parts to derive the last equation and r)n V = diag{l, —1,-1,-1}. 
In (0) G is the gravitational constant. The gravitational field h couples to the gluodynamics by 
means of the second term 

These two terms are invariant under scale and conformal transformations in a space-time of 
arbitrary dimension. However, the scale symmetry in QCD is broken down by non-vanishing 
gluon condensate {F^ U F® U )> 0. This leads to the scale anomaly which manifests itself as an 
anomaly of the QCD energy-momentum tensor (Q^) trace 

0» = ^F« v F^ , (8) 

where the /3-function in the leading order in the strong coupling g reads 



Its vacuum expectation value is 

<0£> = -4|e v | . (10) 

In order for the energy-momentum tensor, calculated from the lagrangian © to satisfy (0), the 
third term in which explicitly breaks the scale invariance has been introduced. This term 
can also be thought of as the cosmological constant. In the absense of the gravitational part (0) 
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of the lagrangian (||) the condition (|i~0| ) is satisfied trivially at the tree level. Performing the 
Legendre transformation of @ one arrives at (|5[) [|| . The mass of the dilaton can be expressed 
through the vacuum energy density and the gravitational constant as follows 

™ 2 = yNG . (11) 

Now we require the energy momentum tensor calculated from the lagrangian (||) to satisfy 
the vacuum normalization condition (|l0|). This means cancelation of the dilaton and gluon 
contributions to the vacuum energy density. We can formulate ( |i~0|) in terms of the vacuum 
expectation values. To this end let us calculate the trace of the energy-momentum tensor 

where the last term in the right hand side is the total derivative § . Equation of motion for the 
dilaton field is 

^ d, (e*/ 2 d, X ) - M (d, X ) 2 + X e* \e v \ - X e*\ F^F^ = . (13) 
Using (^) in ( |I2| ) we arrive at 

9^ = -4\e v \e^ + X e x F^F a ^ . (14) 
Hence, (|To|) is equivalent to the following relation 

( X e x F; v F a ^)=A\e v \(e*-l) . (15) 



It is clear from equations (14) and (13) that in the limit |e v | — ► the classical equation 9£ = 
is recovered. However, quantum corrections can give non-vanishing contributions to (90 in the 
same limit by introducing dimensonal cutoff. Eq. ( |i~5| ) assures that this does not happen in our 
theory. 

Now, we proceed to the derivation of Feynman rules. In addition to usual Fyenman diagrams 
of gluodynamics there are those corresponding to the dilaton propagator, dilaton self-interactions 
and dilaton-gluon interactions. In this work we are interested only in contributions of the order 
of g°. They are shown in Fig. |l|. Other vertices, which contain three and four gluons coupled to 
the dilaton are supressed by factors g 1 and g 2 respectively. In our approach the strong coupling 
constant is necessarily small since all gluon modes which enter the lagrangian (|5|) are hard by 
construction. All soft modes, which carry momenta p < m, have been integrated out. The 
dilaton Feynman propagator can be read from (||): 



m? 



D d (p) = T-r ~ 2 • (16) 

|e v | p z — m z — is 



To completely specify our effective theory we still have to calculate the ultra-violet cutoff Mq. 
Contributions to the vacuum energy density involve loops and therefore they are divergent in 
the ultra-violet regime of the theory. Let us introduce the cut-off Mq which controls divergences 
of vacuum excitations contributing to the vacuum expectation value of the trace of the energy- 
momentum tensor. The same parameter controls divergences in the dilaton phase space. Our 
requirement that the normalization (jl^), or equivalently (|l5|) , holds means that we impose a 
restriction on the possible values of Mq. Consider a contribution stemming from pure dilaton 
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(n - 1) {p\p% - (pi • p 2 )7]^)5 ab 
n = 2, 3, . . . 
(a) 



{( n " !) + ^=T) Ei>j • Pi)} 

n = 3,4, ... 

(b) 



Figure 1: Feynman diagrams of order g° and given number of dilatons n; (a) gluon-dilaton 
interaction, (b) dilaton self interaction. 







(a) 



(b) 



Figure 2: Some vacuum excitations contributing to the energy-momentum tensor trace, (a) 
Dilaton self-interactions; first diagram is proportional to (x 2 ), (b) dilaton-gluon interactions. 



vacuum fluctuations, Fig. |2|(a). In the next section we will argue that for the processes we are 
interested in here the phase space is available for the emission of only two dilatons. In this 
two-dilaton approximation 



(^)vac = -4|e v |(e*> vac « -4|e v | (l + \( X 2 )^ 



Hence, the contribution to the vacuum energy-density is 



2kv|(x 2 



d 3 k 



1 



1 f M o dM 2 m 2 
" 2|£vl Wo 'M 2 M J (27r) 3 2w " 4vr 3 



m 2 M 2 



(17) 



(18) 



Another contribution of the same order in the strong coupling, 0(g ), comes from the first 
diagram shown in Fig. |2](b). We have 



{N 2 c -l)-f 
vr Jo 



M o dM 2 1 



d 3 qi 



M 2 2! J (2tt) 3 2lo 1 J (2tt) 3 2uj 2 



' , "' / - -r 2 - ( Ql ■ q 2 ) V ^f 0(M -m- lo 2 ) 



Ml 



m 



280 • 48 (2vr) ( 



(N£ - 1) — 



(19) 



where T 2 is a two-dilaton phase-space factor, see (|28|) . Contributions (|18| ) and (|19|) cancel out 
when 

/U l\ 2 / 3 

M 2 = Svr 2 2 / 3 105 1 / 3 (N 2 - l)- 1 / 3 ( M ) . (20) 



The pure gluon vacuum diagrams are proportional to the strong coupling constant a s as can 
be easily seen from (g) and (||). Inasmuch as the strong coupling constant is small, those 
contributions to the vacuum energy density can be neglected. 
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3 Derivation of the soft pomeron 



The leading contribution to the total cross section at high energies stems from the ladder dia- 
grams fl9[j . Such a diagram is shown in Fig. ||. Note that the energy dependence of the Born 
amplitude is determined by the spin of particles exchanged in t-channel. Thus the leading contri- 
bution at high energies comes from exchange of vector particles - gluons [20, ^IJ]. These gluons 



prefer to form ladders, since each rung gets enhanced by large logarithm of energy Ins 3> 1. By 
virtue of the optical theorem, the total cross section atot can be calculated as a sum of all such 
diagrams 

^t = EE ff " = -EE Im - M ^ = o)? , (21) 

r>0n>2 S r>0n>2 

i.e. we sum over all possible numbers n > 2 of dilatons produced in s-channel in each rung for 
a given number r of rungs and then sum over all r. 




Figure 3: One of the ladder diagrams contributing to eq. (21 ) 



The most convenient way to calculate the amplitudes associated with ladder diagrams is 
to employ the Weizsacker- Williams approximation. This means replacing almost longitudinally 
polarized Coulomb gluons exchanged in t-channel by transversely polarized (real) ones. In this 
approximation the gluon propagator reads [21]: 



s a q p q qj_ 

where a q and j3 q are Sudakov variables defined by 

q = a q pi + (3 q p 2 + q± ■ 

Weizsacker- Williams approximation is valid as long as the following conditions hold 

1 > a qi > . . . > a qi > . . . » a qr+1 , 
1 » qr+1 > . . . > A« > . . . > (3 qi , 

In the Born approximation (two gluon exchange) there are no produced dilatons. Contribu- 
tion of a one-rung diagram to the total cross section reads 

2 

(27r) 2 8(( Pl - qi ) 2 )8((p 2 +q 2 ) 2 )T n (M) 



(22) 

(23) 

(24) 
(25) 
(26) 



d 4 q 1 f d 4 q 



(2tt) 4 7 (2vr) 



2g 2 ( qi -q 2 ) ± Tr(t a t b 



1 



4/V? 



(27) 
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where T n (M) represents contribution of n dilatons to the phase space and M 2 = (qi — q 2 ) 2 is the 
four-momentum squared available for production in the dilaton vertex. The phase space T(M) 
has been evaluated in Ref. 0: 



r n (M) = 1 (g) n (n - If (2^ 5* (± ki -k)f[ d " h 



J (27T) 3 2Wi 



= ^ J (2?r) (n " 1} n!(n-l)!(n-2)! ' (2§) 

where (n — l) 2 comes from the vertex shown in Fig. |l|(a). Using ( |27| ) and ( |28| ) we derive the 
contribution of one-ladder diagram 

£ ° n r= x = a B -^ir 2 M 2 / ° dM 2 £ r n (M) = ct b Ap In a , (29) 

n>2 V^J JO n=2 

where ctb ~ ct s is a Born amplitude, 

A P = ^ F 2 (-;l,3;a) , (30) 

0-F2 is a generalized hypergeometric function, and the relevant parameter of our theory for high 
energy scattering is 

m 2 Ml , . 

° = M ' <31) 

Expansion of function §F 2 in powers of parameter a 

2 3 

F 2 (-;l,3;a) = 1 + - + — + — + 0{a A ) (32) 
y ' ' ' ' 3 48 2160 y ' v ; 

is equivalent to expansion in the number of emitted dilatons. We will argue below that the series 
(32) are rapidly convergent and actually, there is not enough phase space for emission of large 
number n of dilatons since T n {M) is a rapidly decreasing function of n. Thus, we can safely 
neglect emission of more than two dilatons. 

Now, as we know the contribution of one-rung diagram, calculation of the total cross section 
is straightforward. Indeed, the contributions of different rungs factorize out, leading to the 
Regge-like behavior of the total cross section ( |2l| ) 

00 

(Jtot = cr B ^2 — ln r s A r P = a B s Ap . (33) 
r=0 r ' 

Corrections to the pomeron intercept coming from perturbative QCD are proportional to the 
strong coupling constant g 2 and neglected in this approach. Loop corrections to the dilaton 
propagator give small contributions of the order of 

, „ am 2 

AA *~Mr <<a • 

Finally, upon substitution of © and © into © we obtain the following result for the soft 
pomeron intercept 

ap = l + Ap = l + mfl (N ^ ir ^(^\ i/3 . (34 ) 
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To estimate |e v | in pure gluodynamics we use the fact that it is related to the vacuum energy 
density in QCD |e v |QCD by 

11 . 

M ^ y3|e v | Q cD • (35) 

Indeed, as was argued in Ref. ||, 

(a s F^F a n s iue « 3(a s F^F^) QCD . (36) 

Also the beta function @ in QCD scales with number of col-ours and flavors as /3qcd oc 
(lliV c — 2Nf)/3 whereas in pure gluodynamics (3 g \ ue oc lliV c /3. Since by (||) |e v | oc (3 the 
vacuum energy density in pure gluodynamics gets enhanced by the factor 11/9. Together with 
( |36|) this yields (|35|). Sum rules analysis of || makes it possible to estimate the QCD vacuum 
energy density |e v |QCD = (0.24 GeV) 4 . Owing to ( p5| ) we obtain the following estimate: 



|e v | = 0.012 GeV 4 . (37) 

Another parameter of our effective theory is the dilaton mass. Note, that due to the isospin 
conservation each dilaton can produce only even number of pions in the final state. Let us 
assume that in full QCD with light quarks the dilaton mixes with the broad physical scalar 
<T-resonance [[jj at 

m = 0.4 4-0.5 GeV . (38) 
Note that the spectral density of the dilaton excitations can be spread over several scalar reso- 



nances of different mass (see, e.g., |J3|). However, for our purposes it suffices to keep only the 
lightest resonance, which gives the dominant contribution to the intercept because of the phase 
space constraints. 

Let us now proceed to the numerical estimates. In view of (|37|) and ([^) Eq. ( |20"| ) implies that 
Mq = 3.0 4 2.8 GeV. Substituting these numbers to (|34|) we obtain for the pomeron intercept 



ap = 1.047 4 1.085 which is in reasonable agreement with the phenomenology of high energy 



scattering [22]. Using these numbers in (p2]), we see that the emission of three and more dilatons 
contributes 1 — 0-^2 (— ; 1, 3; 0) /0-P2 ( — ; 1, 3; a)=21 4- 27% to the final result; this can be used as 
an estimate of the accuracy of our two-dilaton approximation. 

We can calculate the inclusive spectrum of dilaton production by observing that Eq. ( p8|) 
is nothing but the number distribution of the dilatons produced per unit of rapidity. Indeed, 
each rung of the ladder diagram corresponds to one unit in rapidity (multiplied by Ap) in the 
leading logarithmic approximation. Define the probability P n (a) to find n dilatons in the final 
state 

*«■)= . (39) 

Then the average number of dilatons produced per unit of rapidity is given by 

= ™ . (40) 
d{Apy) ^ Q F 2 {-;l,3;a) 

From ( |3lT ) we get an estimate of the dilaton phase space parameter a = 0.768 4 1.033. Substitut- 
ing this value into (40) one arrives at (dN x /dy) = 0.1 4- 0.2. Since we have identified the dilaton 



with the physical u-resonance, we have to assume that it decays mostly in two pions. Charged 



8 



pions give approximately 2/3 of the total number of pions. Therefore, the average multiplicity 
of charged pions per unit of rapidity is estimated as 



dN*_ 
{ dy ' 



3 { dy ' 



0.14-^0.26 



(41) 



This agrees with the result obtained in[16|: the soft pomeron is responsible for less then 10% of 
the observed multiplicity. In Fig. || we show the number distribution ( p9| ) of produced charged 
pions. 




n 



pions 



Figure 4: Number distribution of charged pions produced by dilaton mechanism for a = 0.9 at 
hadron-hadron collisions. 



4 Summary and discussion 

Summarizing, we have started with the derivation of the QCD effective lagrangian (|5|) which 
incorporates scale anomaly Q Q and represents propagation and interaction of gluons and dila- 
tons. The dilaton field saturating the sum rules (§) describes the QCD vacuum structure. The 
pure dilaton part of lagrangian (||) (first and second terms) is a low energy effective gluodynam- 
ics obtained by integrating out low energy modes with momenta q such that q < m ~ 0.4 -j- 0.5 
GeV. The third term of (||) represents pure gluodynamics and the interaction of gluon field 
with the dilaton one. It is important as far as momentum transfer does not exceed Mq. It has 
been already pointed out in Ref. || ||, 16] that the scale Mq, determined by the semi-classical 
fluctuations of the gluon fields in vacuum is quite large, Mq ^> A. 

Normalizing the vacuum energy density by (|i"o| ) we obtained an expression (]20| ) for the 
ultra-violet cutoff Mq of the theory in terms of vacuum energy density |e v | and the dilaton 
mass m. Divergent integrals that we have encountered in Eqs. (|l8|),([l9| ) , (l29|) should not surprise 
the reader since the dilaton acts as the gravitational field. Indeed, quantum gravity is a non- 
renormalizable theory since the coupling of gravity to matter fields G has canonical dimension 
(mass) -2 . In technical terms, each emitted dilaton contributes a factor of m 2 /|e v | coming from 
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the normalization of the dilaton propagator which must be compensated for by Mq to maintain 
the scattering amplitude dimensionless. 

In the framework of this effective theory we addressed the problem of soft pomeron. We 
calculated the high energy scattering amplitude in the leading logarithmic approximation by 
summing the ladder diagrams, see Fig. ||. Generally, any number of dilatons can be emitted in 
the i-channel in each ladder rung, but actually most of the phase space is occupied by only two 
of them. Owing to the isospin conservation the dilaton in the final state decays into an even 
number of pions. By identifying the dilaton with the physical a resonance, we conclude that it 
should decay mostly in two pions. We estimate the pomeron intercept ap = 1.0474- 1.085. This 
value agrees well with the analysis of the data on high energy scattering [22]. 

We also estimate the final state multiplicity per unit of rapidity — 0.2. The low 

particle multiplicity resulting from the soft pomeron exchange has been noted already in [16]. It 
implies that the bulk of particles is produced by another mechanism. Such a mechanism must 
involve multiple production of particles per unit of rapidity and perhaps closely related to high 
parton density regime of QCD [21, 24]. 

The reader familiar with the previous work of two of us, Ref. [(Tgl , might be confused by the 
fact that there the opposite to (|34|) dependence of A on the vacuum energy density |e v | and the 
scalar particle (glueball) mass rap has been found. The reason is simple: a glueball mediates 
a short-range strong interaction, so the cross section for its production is proportional to rn^ 4 . 
On the other hand, gravitation is a long-range interaction, the strength of which is proportional 
to masses of interacting particles (dilatons). In other words, the effective lagrangian (||) is 
constructed to obey scale invariance on the classical level. This requirement translates into the 
absence of the coupling of the dilaton to two gluons, because a term linear in the dilaton field 
X in (||) would induce a non-zero vacuum expectation value (x) which would violate the scale 
invariance already on the classical level. (Even though the field x 1S dimensionless, its couplings 
are not, and so its v.e.v. would violate scale invariance of the vacuum state). This is why the 
coupling of the scalar glueball to two gluons, used in [16], does not appear in the approach of 
this paper, which explains the difference in the dependence of the final results on the scalar 
particle mass. 

The reader may also wonder about the meaning of coupling QCD to gravity. One way of 
interpreting this is to say that coupling to gravity is a purely mathematical trick which makes 
it possible to take all symmetries into account in the most elegant way @. 

However, there might be another interpretation which goes beyond the standard model. We 
treated the gravitational constant G as a free parameter. We can find its numerical value using 
(11): G = 0.2 GeV -2 . One might thus assume that gravitational interactions become strong 
already at the scale G -1 / 2 ~2.2 GeV. In this case cutting divergent integrals associated with 
dilaton fields at the scale of Mq means that the discussed theory (|5|) is an effective one. It implies 
that the full renormalizable theory must be recovered at some higher energies. In this sense our 
scale Mq can be derived from some fundamental scale M*. In recent years we have learnt that 
the existence of strong gravity regime at experimentally accessible energies can be realized in 
theoretical models if one assumes that our world has dimension larger then four. The Standard 
Model is confined to the four-dimensional manifold, while gravity propagates in all dimensions. 
One further assumes that there is no specific gravity scale, but rather a unique fundamental 
scale for the Standard Model and gravity around the electroweak scale M* ~ 1 TeV |l^]. Thus 
the observed Plank scale is small since the size of extra dimensions is large (up to mm). This 
does not seem to contradict experimental data since the gravitational potential has been tested 
at distances down to a cm. Existence of strong gravitational interactions at TeV scale has a 
number of remarkable experimental signatures [18. 13, pq |. On the contrary, one may assume that 
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the scale of gravity is different from the Standard Model one, but the gravitational lagrangian 
involves higher derivatives of the metric tensor [27]. In this approach the quantum gravity scale 
can be as low as 10~ 3 eV. However, we explicitly neglected all higher order derivatives of metric 
in (||). Even if the quantum gravity indeed becomes strong at M* <C Mpi the probability that 
a dilaton escapes into higher dimensions is suppressed by factor of order of 0(s/M* 2 ), which 
means that our four-dimensional effective lagrangian (||) is a good approximation at energies 

<C M*. So, at these energies gravity can be considered as a background classical field, 
since quantum corrections are small. It would be interesting to address the question of what 
is the mechanism by which the full theory at the fundamental scale yields the conformally flat 
GeV-scale gravitational effects. However, we leave this problem aside. 

Whatever the underlying mechanism is, the broken scale invariance can be encoded in the 
effective low-energy theory of gluodynamics in the form of the lagrangian d). We hope that it 
provides a useful tool for a systematic approach to non-perturbative processes at high energies. 
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